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Abstract 


What can we learn from solar neutrino observations? Is there any solution to the 
solar neutrino anomaly which is favored by the present experimental panorama? 
After SNO results, is it possible to affirm that neutrinos have mass? In order to 
answer such questions we analyze the current available data from the solar neu¬ 
trino experiments, including the recent SNO result, in view of many acceptable 
solutions to the solar neutrino problem based on different conversion mecha¬ 
nisms, for the first time, using the same statistical procedure. This allows us 
to do a direct comparison of the goodness of the fit among different solutions, 
from which we can discuss and conclude on the current status of each proposed 
dynamical mechanism. These solutions are based on different assumptions: (a) 
neutrino mass and mixing, (b) non-vanishing neutrino magnetic moment, (c) 
the existence of non-standard flavor-changing and non-universal neutrino inter¬ 
actions and (d) the tiny violation of the equivalence principle. We investigate 
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the quality of the fit provided by each one of these solutions not only to the total 
rate measured by all the solar neutrino experiments but also to the recoil electron 
energy spectrum measured at different zenith angles by the Super-Kamiokande 
collaboration. We conclude that several non-standard neutrino flavor conversion 
mechanisms provide a very good fit to the experimental data which is compa¬ 
rable with (or even slightly better than) the most famous solution to the solar 
neutrino anomaly based on the neutrino oscillation induced by mass. 

PACS numbers: 14.60.Pq, 13.15.+g, 26.65.+t 


Typeset using REVTgX 


2 



I. INTRODUCTION 


Solar neutrino observations coming from Homestake |T|, Karniokande J2]] r SAGE [.^], 
GALLEX [|], GNO and Super-Karniokande (SK) [6] have been suggesting a picture which 
is conflicting with the predictions from the Standard Solar Model (SSM) (7[-|9|], strongly indi¬ 
cating disappearance of solar electron neutrinos on their way from the Sun to the terrestrial 
detectors. This has been known for many years as the Solar Neutrino Problem (SNP) [7J. 

The extraordinary new result from the Sudbury Neutrino Observatory (SNO) m , inau¬ 
gurates a new era in the quest for the solution to the long-standing puzzle of missing solar 
neutrinos. For the first time in the history of solar neutrino observations, a direct indication of 
the presence of non-electron active neutrino component in the solar neutrino flux is obtained. 
This cannot be explained by any conceivable modification of the SSM but do require some de¬ 
parture from the standard electroweak theory. The indication of non-electron active component 
is based on the difference of 8 B neutrino flux detected through charged current events in SNO 
and the neutrino electron elastic scattering events observed by the SK collaboration, the former 
obtaining a lower rate than the latter. Such difference can be explained by the conversion of 
electron neutrinos into active non-electron (z/ M or u T ) neutrinos along their trajectory from the 


Sun to the detectors at the Earth 10,11 


As we will see, the hypothesis of no flavor conversion of solar electron neutrinos is strongly 
in conflict with the prediction of the SSM; it is now only acceptable at very small confidence 
level (~ 10 -12 ) if only the total rates from solar neutrino experiments is considered (see Sec. 
IV). The disagreement among the observed solar neutrino data and the theoretical predictions 
can be relaxed to 4 a level (7 x 10 -5 ) |12| if one allows all the solar neutrino fluxes to be 
free parameters in fitting the measured solar neutrino event rates. However, this can only be 
obtained under the extreme assumption of vanishing 'Be neutrino flux, which is quite difficult 
to explain. 

Several mechanisms can induce neutrino flavor conversion when one assumes that neutrinos 
are endowed with some properties not present in the minimal standard electroweak theory |H|] . 
The most well known mechanism is the neutrino oscillation induced by mass and mixing JT^-|T^]. 
The fact that the terrestrial experiments are less sensitive to these resulting non-electron neu¬ 
trinos can explain their observed lower counting rates. The purpose of this article is to compare 
quantitatively the capabilities of several different mechanisms to explain solar neutrino data, 
for the first time, based on the same statistical procedure. 

Readers are invited to take a look at Tables m and 13 which summarize our most important 
results. A number of possible solutions to the solar neutrino anomaly still survives even after 
SNO results, fitting the data with significant confidence level. In fact, combined analysis of 
the data, which includes not only the total rates measured by all solar neutrino experiments 
but also some information which is independent of the total neutrino flux, namely, the energy 
spectrum and the zenithal dependence of the data, suggests that the large mixing angle MSW 
solution in matter, as well as the mechanisms based on resonant spin-flavor precession, non¬ 
standard neutrino interactions and violation of the equivalence principle all provide a fit of the 
data with the conhdence level > 60%. 
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In Sec. [H. we briefly review several mechanisms that induce flavor conversion of solar neutri¬ 
nos, which will be discussed in this work: (a) mass induced oscillation in vacuum and in matter 
(b) resonant spin-flavor precessions induced by a non-vanishing neutrino magnetic moment, (c) 
the existence of non-standard neutrino interactions inducing flavor-changing and non-universal 
currents and (d) the violation of the equivalence principle. In Sec. B the procedure of the 
statistical analysis used in this work is presented, while our results are given in Sec. B Finally, 
Sec. 0 is devoted to the discussion of our results. 


II. NEUTRINO CONVERSION MECHANISMS 


The neutrino conversion mechanisms we will consider in this work can be phenomenologi¬ 
cally described by the following Schroedinger-likc evolution equation: 
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where r is the distance travelled by neutrino, v e is the initial electron neutrino state, v is the 
neutrino state to which the conversion mechanism leads and the explicit form of the elements 
of the 2x2 Hamiltonian matrix, A{r ) and C'(r), which are in general depend on the position 
r, will be given in the next subsections for each conversion mechanism. For simplicity, we 
assume neutrino conversion only between two flavors, or more precisely, we assume that neutrino 
conversion relevant for solar neutrinos can be effectively described, in good approximation, in 
terms of two flavor conversion. 

In all the mechanisms, either in vacuum or in matter, the phenomenon of neutrino mixing 
occurs and this can be generally expressed as 

V e 
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cos 9 sin 9 
— sin 9 cos 6 
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( 2 ) 


where 6 is the mixing angle which relate neutrino flavor and the propagating eigenstates v\ and 
v 2 hr matter or in vacuum. Note that 9 can be defined even for massless neutrinos as in the 
case of NSNI or VEP. 

For all the mechanisms we analyze in this work, we do not consider the case where v is a 
sterile (or electroweak singlet) neutrino [18] because of the present indication of the presence 
of non-electron active component in the solar neutrino flux, which is provided by the combined 
information from SNO and SK results. Namely, neutrino conversion of solar v e only into a 
sterile neutrino is not favored by the current data. 


A. Mass induced oscillation (MIO) 


Let us first consider the most popular conversion mechanism, the one induced by neutrino 
mass and mixing. If one assumes that neutrinos are massive, the flavor eigenstates do not 
coincide, in general, with the mass eigenstates, leading to neutrino flavor oscillation in vac¬ 
uum [ 14,151 as well as in matter where it can be resonantly enhanced [ 16,17|. The latter is 
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known as the Mikheyev-Smirnov-Wolfenstein (MSW) effect. It is generally believed that this 
is the most plausible mechanism which can induce neutrino oscillation since to introduce mass 
and mixing in the leptonic sector is the simplest and most natural extensions to the standard 
electroweak model, so theoretically well motivated. 

In both scenarios, MIO can convert solar electron neutrinos into neutrinos of a different flavor 
and consequently explain the deficit of observed solar neutrino with respect to the predictions 
of the SSM. We will refer to these type of solutions as standard solutions to the SNP, they are 
very well described in recent Refs. [19,20]. MIO requires in Eq. ( 0 ) that 
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C(r) 


5m 2 G F . 

■—cos29 +75 iV e (r), 
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( 3 ) 


sin 26, 


where dm 2 = — m\ is the mass squared difference of the two neutrinos involved, 6 is the 

vacuum mixing angle, G F is the Fermi constant, E the neutrino energy and N e {r) is the electron 
number density at position r, and here v is identified with or u T (or their linear combination). 
We compute the conversion probability using the analytic formulas, properly taking into account 
the neutrino production distributions, as well as the Earth matter effect as in Ref. |H[]. The 


relevant oscillation parameters, which must be determined by the fit to experimental data, are 
5m 2 and 9. Recent post-SNO analyses can be found in Ref. B Depending on the allowed 
parameter values, the solutions based on this mechanism are classified into large mixing angle 
(LMA) MSW, small mixing angle (SMA) MSW, low-<5m 2 (LOW) MSW and vacuum oscillation 
(VAC) solutions. 


B. Resonant spin-flavor precession (RSFP) 


Assuming neutrinos have a non-vanishing transition magnetic moment, electron neutrinos 
interacting with the solar magnetic field can be spin-flavor converted into active non-electron 


anti-neutrinos [[L|, if they are of Majorana type. Here we do not consider the case of Dirac type 
neutrinos since it involves a sterile one. Furthermore, such spin-flavor precession of neutrinos 
can be resonantly enhanced in matter [EM|, in close analogy to the MSW effect [|TC] . [TT|| . RSFP 
could strongly depend on the neutrino energy and provoke different suppressions for each por¬ 
tion of the solar neutrino energy spectrum, in a similar way as in the case of the MSW effect. 
In fact it has been known that RSFP could provides a satisfactory description p5|-[27[| of the 
actual experimental panorama [|I|-|5[] . However, see a recent work [|2H| where it is found that the 
non-resonant spin-flavor precession can also explain the solar neutrino data. 

The time evolution of neutrinos interacting with a magnetic held B through a non vanishing 
neutrino magnetic moment /i v in matter is governed by Eq. (0) by defining A{r) and C(r ) as 

m- 


Sin 2 1 

A(r) = -— + ^=G F [N e (r)-N n (r)}, 
C(r ) = fJ. v B(r), 


( 4 ) 
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where N n (r) is the neutron number densities, /j, u is the transition magnetic moment between two 
neutrino involved, B(r ) is the solar magnetic held and here v is identified with or U T (or their 
combination). We assume, for simplicity that the mixing angle is zero in this mechanism. Note 
that, roughly speaking, fi u B(r) is playing the same role as the mixing term (<5m 2 /4i7) sin 29, 
which appeared in the standard mass induced oscillation mechanism. In this work, we assume, 
as a reference value, the magnitude of the neutrino magnetic moment to be /i„ = lO -11 ^ (//# 
is the Bohr magneton). Since the relevant quantity is only the product of the neutrino magnetic 
moment and the magnetic held, it must be understood that if /i„ is taken to be smaller, the 
solar magnetic held must be properly increased to achieve the same effect. 

The RSFP mechanism crucially depends on the solar magnetic held prohle along the neu¬ 
trino trajectory. In our present analysis we assume a particular typpe of solar magnetic held 
prohle which could explain well the solar neutrino data (before SNO result) with relatively 
weak magnetic held (~ 0(10) kG), which was considered in number of works in Refs. mm- 
The prohle we will use has a triangular configuration in the solar convective zone. For definite¬ 
ness, we take the one which is found in Ref. [^{J. The prohle has a vanishing magnetic held 
in the internal part of the Sun, linearly growing from r = 0.7 to r = 0.85 where it achieves its 
maximum value R max and begins to linearly decline until the surface of the Sun, where r — 1. 
Here r is the radial distance from the center of the Sun normalized by the solar radius. In this 
work, we consider B max up to 500 kG taking into consideration the upper limit (300 kG) of the 
magnetic held at the bottom of the convective zone obtained in Ref. p0|. Once the shape of 


the magnetic held prohle is hxed, the relevant parameters which must be determined by the ht 
can be chosen to be 5m 2 and B max . 

We compute the conversion probability by numerically integrating the evolution equation 
as in Ref. |27|. For simplicity, we assumed that all the neutrinos are created in the solar center 


and neglected the Earth matter effect, as discussed in Ref. p7 


C. Non-standard neutrino interactions (NSNI) 


In his seminal paper Wolfcnstein ||T6[ observed that NSNI with matter can also generate 
neutrino oscillation even without havor mixing in vacuum. Some explicit examples of such 
NSNI induced neutrino oscillation were considered in Refs. ||RT| . pT2| . It has been shown that 
NSNI could be relevant to solar neutrinos propagating in the solar matter along their path 
from the core of the Sun to its surface |3T H33l1 as well as in the Earth until they reach the 
detector. 

The evolution equations for massless neutrinos (or neutrinos with degenerate mass) hav¬ 
ing such NSNI in matter can be phenomenologically expressed by Eq. ([!]) with the following 
dehnitions of A(r ) and C{r ) 
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A{r) = T GF [JV e (r) - 
C(r) = V2 G F e,N,(r), 
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where e'f and €f are the phenomenological parameters which characterize the strength of the 
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NSNI with fermion / whose number density is given by Nf(r) with / = u or d quark, and here 
v is identified with or v T (or their combination). 

Here, by taking into account the charge neutrality, Nf can be written in terms of electron 
and neutron number densities, as follows, 


JV/( r) 


JV„(r) + 2/V„(r) (/ = «) 
2JV„(r) + N e (r) (f = d). 


( 6 ) 


The parameters 6f appearing in the off-diagonal elements of the Hamiltonian matrix is re¬ 
sponsible for flavor changing neutrino interactions, it plays a similar role to the mixing term 
(5m 2 /AE) sin 29 in the MIO mechanism. It exists even if there is no neutrino mixing in vacuum. 
On the other hand, e'f appearing in the diagonal element which is responsible for flavor diag¬ 
onal neutrino interactions with matter |JTJ, somehow plays a role similar to the mass squared 
difference in the MIO mechanism, since this term leads to resonantly enhanced conversion when 
its magnitude coincides with that of the standard electroweak neutrino interactions at some 
point r rea along the neutrino trajectory, satisfying the resonance condition, 


s'fN f (r ies ) = N e (r res ). 


(7) 


An immediate consequence is that if / was to be identified with electron, it cannot induce 
resonant neutrino flavor conversion |34[ . 


An important characteristic of this mechanism is that the conversion probability does not 
depend on neutrino energy as it is understood from the evolution equation. Nevertheless this 
mechanism can explain quite well the solar neutrino data which imply strong energy dependent 
conversion. The reason is that even though the conversion probability itself is completely energy 
independent, after taking average of the probability over production distributions, different 
neutrinos from different nuclear reactions in the Sun (pp, 8 B, 'Be, etc.) can have different final 
average probability due to the fact that their production distributions are different, provided 
that resonant conversion occur close to the solar core (~ 10% or so of the solar radius). 

Following Ref. |33|], we first compute the conversion probability using the analytic formulas 
for a given production point in the Sun and then take average over the production distribution 
for each neutrino source. We also take into account the Earth regeneration effect, which is 
important for some region of the parameter space. 

As it was observed in Ref. |JT|, only v e —> v T conversion are compatible with the existing 
phenomenological constraints on e/ and e'f. In this way, the relevant parameters which must 
be determined by the fit for this mechanism are e'f and e/. 


D. Violation of the equivalence principle (VEP) 

It has been proposed that violation of the equivalence principle could induce neutrino flavor 
oscillation even if neutrinos are massless (35|.|3(J. In this scenario, neutrino mixing and flavor 
oscillation can be induced if two (or more) neutrinos involved have different gravitational cou¬ 
plings which imply VEP. In this case, weak interacting eigenstates and gravitational interacting 
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eigenstates can be different and will be related by a unitary transformation that can be pa¬ 
rameterized, assuming only two neutrino flavors, by a single parameter, the mixing angle 0 G 
similar to the case of neutrino mixing in vacuum induced by mass. 

The evolution equations for these flavors, which are assumed to be degenerate in mass, 
propagating through the gravitational potential 0(r) in the absence of matter is given by Eq. (1|) 
with p5| , p6| : 


A(r ) = 2E(p(r)5'y cos 29 g , 
C(r) = 2E(j)(r)5'ysm29 G , 


( 8 ) 


where is the quantity which measures the magnitude of the violation of the equivalence princi¬ 
ple, the difference of the gravitational couplings between the two neutrinos involved normalized 
by the sum. ffere 2E<p{r)8^s\n26 G plays a similar role to the mixing term (8m 2 /AE) sin 26 in 
the MIO mechanism, this means the VEP mechanism gives rise to an oscillation length inversely 
proportional to the neutrino energy E, while MIO expects it to be directly proportional to E. 

As in Ref. |57|] , we take 0(r) to be constant (~ 10 -5 ), assuming that the local supercluster 
contribution is the dominant one |38|. The relevant parameters which must be determined 
by the fit in this mechanism are |05y| and 0 G . We consider the product of 0 and since 
the former has a large uncertainty and only the product is relevant in the fit. Similar to the 
mass induced oscillation, we can have two types for the VEP mechanism: (a) VEP induced 
MSW-likc resonant conversion and (b) VEP vacuum conversion. The former was first discussed 
in Ref. [36], and then analyzed by several authors p9 


However, Ref. 


showed that it is 


not a good solution because the required magnitude of the parameters is incompatible with 
the CCFR experiment results |TT|] which exclude \(j>8 , y\ larger than ~ 10 -23 . For the latter one, 
it has been shown that this mechanism is a new solution for the solar neutrino anomaly [[37 
Recent analysis of VEP vacuum solution was discussed in Ref. m- 


We calculate the conversion probability using the analytic formulas as in Ref. |37|]. Sim¬ 
ilar to the case of vacuum oscillation solution in the MIO scenario, we can neglect, as good 
approximations, the neutrino production distributions and the Earth matter effect [[I7j . 


III. STATISTICAL PROCEDURE 


Our main goal is to determine the allowed values of the relevant parameters for each one 
of the mechanisms that can explain the experimental observations without modifying the SSM 
predictions and to evaluate the quality of the fit. In order to achieve our goal, we adopt the 
minimum y 2 statistical treatment of the data following the description found in Ref. [|43| which 
was also employed in Refs. | PT| . pT7| . ^TT] . po| . pT| . [T~F|| but with some modifications. 

In this work, we perform three kinds of analysis using three y 2 functions, Xr f° r the analysis 
of the rates only, y| for the flux independent analysis, and y 2 omb = y 2 R + y| for the combined 
analysis, where the definitions of y 2 R and y| are given in the next subsections. For each case, 
we minimize the y 2 function in order to determine the best fitted parameters as well as its 
minimum value yV n , which is relevant for the evaluation of the goodness of fit. The allowed 
parameter region can be determined by the condition, 


8 












































( 9 ) 


X 2 < Xmin + A \' 2 , 

where Ay 2 = 4.61, 5.99 and 9.21 for 90%, 95% and 99% C.L., respectively. 


A. y 2 for the analysis of the rates 


First we describe the y 2 for the analysis of the total event rate measured by the Chlorine 
(Cl) experiment |l|, the Gallium detectors GALLEX/GNO f||| (we use the combined results 
of GALLEX and GNO) and SAGE [j3|, the water Cherenkov experiment SK [p|] and also the 
recent neutrino-deuterium charged current data by SNO (10|. For simplicity we do not consider 
the result from the Kamiokande experiment || as it is consistent with the current SK one and 
has much larger experimental errors. All the solar neutrino data for the analysis of the rates 
used in this work are summarized in Table |I|. In all our calculation of the theoretical predictions 
we have used the latest “BP2000 SSM+New 8 B” p|| fluxes, which include the recent the new 
measurement of Si 7 (0) [|IS|. Our y 2 function is defined as follows: 


xk = E 


E>th 7> 

K l ~ K , 


obs 


i -1 r 
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pth pi 
K J ~ R J 


obs 


( 10 ) 


where R\ h and R° hs denote, respectively, the predicted and the measured value for the event 
rates of the Eve solar experiments considered, i = Cl, GALLEX/GNO, SAGE, SK, SNO. 

In order to compute the predictions for the rates, for Chlorine, Gallium and SK detectors, 
we follow Refs. For SNO, we compute the rates using the neutrino charged current 

cross section on deuterium, u e + d —> p + p + e~, as given in Ref. [fhj and taking into account 
the energy resolution as described in Ref. [ICj. The error matrix ctr contains both experimental 
(systematic and statistical) and theoretical errors |H3| . 


B. y 2 for flux independent analysis: SK zenith and spectrum 


Our final result is the fit derived from the combined analysis of all presently available solar 
neutrino data, which includes the flux independent information presented by SK Collaboration 
47|. We include this information using the SK data split simultaneously on 7 zenith bins 


and 8 spectrum bins. We will refer to these as the SK zenith-spectrum data. These data 
are summarized in Table |TI|. For the first and last spectrum bins the zenith data are grouped 
together in one data point, due to low detection rate at these bins, in total 44 data points related 
to flux independent information. We took the statistical and systematic energy-correlated and 
uncorrelated uncertainties, listed in Ref. [^]. The systematic uncertainties are assumed to be 
fully correlated in the zenith angle splitting. We then use the following y 2 expression for the 
flux-independent information: 


Xl = 


— -^SK,J [^(h j) 2 ] 1 [ a -RsK,i — R SK,i\ ) 


->obs 


->th 


■>obs 


(ii) 


i=l,44 


where Rg 1 ^ i is the theoretically expected event rates for the i-th bin computed by using the 8 B 


neutrino energy spectrum given in Ref. B3] normalized to the BP2000 SSM value 


^SK,i B 
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the corresponding observed rate reported by the SK Collaboration ||7]], a is a free parameter 
to avoid double-counting of the SK total rate in the statistical treatment, and is the 44 

x 44 error matrix for SK zenith-spectrum data. This treatment is slightly different from the 
one presented in recent solar neutrino analysis |9l. |TT| . |22| which used the SK elastic scattering 
spectrum data taken during the day and the night separately whereas we used the SK spectrum 
measured at different zenith angle (night spectrum data are divided into 6 zenith bins). 


IV. RESULTS 


In this section, we discuss the results of our statistical analysis. In Table |T|, a comparison of 
the existing solutions to the SNP, when only the total rates are taken into consideration in the 
statistical analysis is presented. For each one of the indicated mechanisms the best fit values of 
the relevant parameters are shown, followed by the corresponding Xmm and its confidence level. 
In the first row one can find the result of the poor fit of the SSM predictions to the data : xLm 
= 62.5 for 5 d.o.f. which represents a confidence level of only ~ 10 -12 . All the other solutions 
present 2 free parameters to fit 5 experimental data points, resulting in 3 d.o.f.. The largest 
confidence level (87%), showed in Table 0 of the fit to the observed data is achieved by the 
solutions based on neutrino NSNI and the second largest one (80%) is achieved by the RSFP 
solution. We note that the RSFP solution has an extra freedom in choosing the solar magnetic 
held profile. In deed we found that the RSFP solution is rather sensitive to the relatively small 
change of the magnetic held prohle. We found that if we use the similar but different magnetic 
held prohle used in Ref. |27|], which has the peak at r = 0.8, somewhat inner region of the sun, 
solutions exist only for R max larger than ~ 100 kG. Among the solutions with no such extra 
freedom, MIO in vacuum is the second best one (45 %). 

For the no-oscillation hypotheses, in the case where only the SK zenith-spectrum information 
is used, we obtain = 37.7 with a = 0.446, for 43 d.o.f., which is compatible with the 
experimental data at 70% C.L.. This results imply that there are no strong distortions in the 
the shape of SK zenith-spectra when compared to the SSM predictions, and the observed data 
are consistent with no-oscillation hypothesis provided that the overall rate is normalized. 

When we include in the statistical analysis not only the total rates but also the energy spec¬ 
trum and zenithal dependence of the data coming from the SK experiment, then the panorama 
changes as summarized in Table 0. As in Table 0. the first row here shows the fit of the SSM 
predictions to the combined data: xLm = 100.0 for 48 d.o.f., which represents a confidence level 
of ~ 10~ 5 . The solutions listed in Table 0 present 2 free parameters to fit 48 experimental 
data points, resulting in 46 d.o.f.. Here, the large mixing angle MSW solution based on MIO, as 
well as RSFP and NSNI solutions are equally probable all having > 80% C.L.. YEP solution 


also provides a very good fit, 60 % C.L. 

In Figs. l(a)-5(a), we present, for all the conversion mechanisms, (i) the allowed parameter 
region determined from the total rates only and (ii) the excluded parameter region determined 
only from the SK zenith-spectrum information. In Figs. l(b)-5(b), we present the allowed 
parameter region determined from the combined data of the rates and the SK zenith-spectrum 
information. 
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It is worthwhile to note that for the MIO and NSNI cases, Figs.[l](a), 0(a) and |J(a), in some 
parameter regions, there is a significant overlap between the region allowed by rates only and 
the one excluded by the SK zenith-spectrum information. Namely, there is a strong conflict 
between the rates and the SK zenith-spectrum fit for some parameters and this is the reason 
why the allowed region decreases significantly when both data are combined as we can see in 
Figs. fj(b), |](b) and 0(b). On the other hand, for the RSFP and VEP cases, there is no such 
kind of strong conflict between rates and and SK zenith-spectrum fit and therefore, the allowed 
region for the combined data is rather similar to the one determined only by the rates as we 
can see in Figs. 0 and 0. 


V. DISCUSSIONS 


Even though the conventional mass induced oscillation mechanism, which is theoretically 
well motivated, can be considered as the most plausible solution to the solar neutrino problem, it 
is important to realize that solutions based on New Physics in the neutrino sector, such as large 
neutrino magnetic moment, neutrino flavor-changing and non-universal processes, violation of 
the equivalence principle, can still be viable ones providing a fit to the solar neutrino data 
which is comparable to the one based on the conventional mass-induced neutrino oscillation, 
as we showed in this work. 

We also note that some of these mechanisms which do not require neutrino mass in the 
fit to the solar neutrino data do have some close relation with neutrino mass generation. For 
example, in our phenomenological approach, we ignored neutrino mass in the solution based 
on non-standard neutrino interactions but there is no available model that prevents neutrinos 
to acquire mass at radiative level, if flavor changing and flavor non-universal interactions with 
quarks are present, even if no tree-level mass terms appear in the model. 

It is fundamental to mention that it is difficult to explain the atmospheric neutrino prob¬ 
lem [[Ej as well as the LSND anomalies [Q by these alternative mechanisms [[H]]. Furthermore, 
while mass induced oscillation and resonant spin-flavor precession solutions to the SNP can be 
easily conciliated with the standard neutrino oscillation solution to the atmospheric neutrino 
problem, it is not a trivial task to answer if the non-standard flavor-changing and non-universal 
neutrino interactions and violation of equivalence principle solutions to the SNP modify or even 
damage this standard solution to the atmospheric neutrino anomaly |52| . 


Having this picture in mind we conclude that the no specific solution is preferred by the 
current solar neutrino data, although some solutions may have difficulties in reconciling atmo¬ 
spheric neutrino observations. We emphasize that the solar neutrino observations alone can 
not yet conclude if neutrinos have non-vanishing mass or magnetic moment. Furthermore, no 
stringent limit on the existence of NSNI nor on VEP can be currently set by the solar neutrino 
data. 

The ultimate goal of the solar neutrino observations is, of course, to perform a direct exper¬ 
imental identification of the solution. For this purpose, we will have to wait for the upcoming 


11 









solar neutrino experiments, which hopefully will provide a lot of new informations, to reveal 
the true nature of the flavor conversion mechanism which is behind the SNP. For instance, low 
energy solar neutrino experiments, such as Borexino 0, can possibly discriminate among the 


0 


55 


may measure re¬ 


solutions considered in this work, as can be seen in Fig. 7 of Ref. 

In the near future, a new reactor experiment, KamLAND 
disappearance. For all mechanisms studied here, with exception of LMA, we expect no sig¬ 
nificant disappearance in KamLAND because the baseline is too short to developed oscillation. 
Then a positive evidence for neutrino oscillation in KamLAND can establish LMA solution. 
For the negative evidence no conclusion can be drawn to favor a particular mechanism studied 
in this work. A discussion about the consequences of KamLAND experiment on the mechanism 
of neutrino mass and mixing hypothesis can be found in Ref. |[A] . 


Summarizing our conclusions, 

1. The present solar neutrino data only by themselves can not discard, a priori, any of 
the solutions discussed here: (a) neutrino mass and mixing, (b) non-vanishing neutrino 
magnetic moment, (c) the existence of non-standard flavor-changing and non-universal 
neutrino interactions and (d) the tiny violation of the equivalence principle. We refer 
to Tables H and E3 for comparison. All solutions have confidence level over 60% C.L. 
providing a very good fit to the solar neutrino data; 

2. A very robust statement is that LMA is the best preferred solution for mass induced 
oscillation scenario, whereas the SMA is the worst one; 

3. Future experiments could test these different scenarios and possibly discard some of them. 
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TABLES 


TABLE I. Observed solar neutrino data used in this analysis presented together with the theoret¬ 
ical predictions of “BP2000 SSM+New 8 B” g||. 


Experiments 

Observed Rates 

SSM Predictions 

Units 

Homestake 

2.56 ± 0.23 

8-59lii 

SNU 

SAGE 

77.0 ± 6. ± 3 

1301? 

SNU 

GALLEX + GNO 

73.9 ± 4.7 ± 4.0 

1301? 

SNU 

SK 

0.391 ± 0.014 

1 no +0 - 14 

1 - uu -0.15 

5.93 x 10 6 cm~ 2 s _1 

SNO CC 

0.296 ± 0.024 

1 00 +0 ' 14 
i - uu -0.15 

5.93 x 10 6 cm 2 s 1 


TABLE II. Observed SK zenith-spectrum event rates in unit of the BP2000 g, taken from Table I 
of Ref. [47]. Errors are only statistical. For the energy-uncorrelated and energy-correlated systematic 
uncertainties, see Table I in Ref [47]. 


Observed rates and statistical errors in units of SSM 
Day Mantle 1 Mantle 2 Mantle 3 Mantle 4 Mantle 5 Core 


cos 6 Z 


-1.00-0.00 0.00-0.16 0.16-0.33 0.33-0.50 0.50-0.67 0.67-0.84 0.84-1.00 


E e (MeV) 
5.0-5.5 


5.5- 6.5 

6.5- 8.0 
8.0-9.5 

9.5- 11.5 

11.5- 13.5 

13.5- 16.0 


0.431T0.022 0.464T0.060 
0.461T0.013 0.524T0.036 
0.437T0.014 0.449T0.038 
0.434T0.015 0.432T0.042 
0.456T0.026 0.496T0.071 
0.482T0.056 0.532T0.155 


16.0-20.0 


0.436T0.046 
0.410T0.055 0.442T0.048 
0.506T0.033 0.438T0.028 
0.482T0.036 0.460T0.031 
0.493T0.040 0.446T0.034 
0.290T0.055 0.394T0.053 
0.775T0.171 0.685T0.141 
0.476T0.149 


0.453T0.048 

0.466T0.027 

0.503T0.031 

0.448T0.034 

0.477T0.056 

0.607T0.130 


0.495T0.054 

0.424T0.030 

0.461T0.034 

0.435T0.037 

0.439T0.061 

0.471T0.128 


0.434T0.058 

0.409T0.033 

0.439T0.037 

0.484T0.044 

0.465T0.068 

0.539T0.153 
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TABLE III. Comparison of the existing solutions to the SNP when only the total rates are taken 
into consideration in the statistical analysis. For each one of the indicated mechanisms the best fit 
values of the relevant parameters are shown, followed by the corresponding Xmin an< i its confidence 
level. In the first row we see the result of the poor fit of the SSM predictions to the data. All other 
solutions present 2 free parameters to fit 5 experimental data points, resulting in 3 d.o.f.. 


Mechanism 



v 2 . 

Wrnin 

C.L. 

SSM with no oscillation 



62.5 (5 d.o.f) 

4 x 10 -12 

MIO 

dm 2 (eV 2 ) 

tan 2 6 



Vacuum 

6.80 x 10 -11 

0.425 

3.02 

39% 

SMA 

7.91 x 10 -6 

1.59 x 10 -3 

4.60 

20% 

LMA 

2.80 x 10 -5 

0.320 

2.62 

45% 

LOW 

1.05 x 10 -7 

0.743 

10.5 

1.5% 

RSFP 

8m 2 (eV 2 ) 

-Bmax(kG) 




1.11 x 10 -8 

338 

1.01 

80% 

NSNI 

e' 

e 



d-quarks 

0.599 

3.22 x 10 -3 

0.72 

87% 

u-quarks 

0.428 

1.39 x 10" 3 

0.73 

87% 

YEP 

I^Aql 

sin 2 2 0 G 




1.56 x 10 -24 

1.0 

6.23 

10% 
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TABLE IV. Same as in Table III but for the statistical analysis taking into consideration not only 
the total rates but also the energy spectrum and zenithal dependence of the data. 


Mechanism 



Xmin 

C.L. 

SSM with no oscillation 


100.0 (48 d.o.f) 

1.6 x 10 -5 

MIO 

dm 2 (eV 2 ) 

tan 2 0 



Vacuum 

4.65 x 10” 10 

1.89 

46.1 

47% 

SMA 

4.93 x 10" 6 

4.35 x 10“ 4 

61.5 

6.3% 

LMA 

6.15 x 10" 5 

0.349 

38.7 

75% 

LOW 

1.01 x 10" 7 

0.783 

45.0 

38% 

RSFP 

<5m 2 (eV 2 ) 

£W(kG) 




1.22 x 10 -8 

440 

38.4 

78% 

NSNI 

e' 

e 



d-quarks 

0.599 

3.23 x 10" 3 

37.9 

80% 

u-quarks 

0.428 

1.40 x 10" 3 

37.9 

80% 

VEP 

I0A 7 | 

1.59 x 10~ 24 

sin 2 2 6q 

1.0 

42.9 

60% 
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FIG. 1. In (a) the allowed region for the MIO solutions to the SNP at 90, 95 and 99 % C.L. with 
rates only as well as the excluded region at 99 % C.L. from the SK zenith-spectrum information are 
shown. In (b) the allowed region from the combined information from rates and zenith-spectrum data 
is presented. 
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Excluded by SK Zenith-Spectrum at 99 % C.L. 


(a) Rates Only vs. SK Zenith-Spectrum 
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FIG. 2. Same as Fig. || but for the RSFP solution in the B max — dm 2 plane. 
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(b) Rates + SK Zenith-Spectrum 





FIG. 3. Same as Fig. |l] but for the solution based on NSNI with d-quarks in the parameter space 
of e(= e d ) and e'(= e' d ). 
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FIG. 4. 

e'(= e' u ). 


Same as in Fig. || but for NSNI with u-quarks in the parameter space of e(= e u ) and 
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FIG. 5. Same as Fig. |l] but for the VEP solution in the sin 2 2 6q — |0<Fy| plane. Notice in the upper 
panel the very small exclusion region at 99 % C.L. from the SK zenith-spectrum information appeares 
at sin 2 29g ~ 1 and \(/)S^\ ~ 10~ 24 . 
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